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The heavy-fermion metal YbRti2Si2 is studied by low-temperature magnetization M(T) and specific-heat 
C(T) measurements at magnetic fields close to the quantum critical point (H c = 0.06 T, H _L c). Upon ap- 
proaching the instability, dM/dT is more singular than C(T), leading to a divergence of the magnetic Griineisen 



ratio T, 



-(dM/dT)/C. Within the Fermi liquid regime, r mag = -G r (H - h{") with G, = -0.30 ± 0.01 
and H J C " = (0.065 ± 0.005) T which is consistent with scaling behavior of the specific-heat coefficient in 
YbRh 2 (Sio.95Ge .o5)2- The field-dependence of dM/dT indicates an inflection point of the entropy as a function 
of magnetic field upon passing the line T*(H) previously observed in Hall- and thermodynamic measurements. 



mag 



PACS numbers: 71.10.HF,71.27.+a 

Quantum criticality is a topic of extensive current research 
interest in condensed matter physics because it results in un- 
usual finite-temperature properties and promotes the forma- 
tion of novel states like unconventional superconductivity (for 
recent reviews see |l|,|2|,|3j]). A quantum phase transition oc- 
curs at T — upon tuning an external parameter r like pres- 
sure, doping or magnetic field to a critical value r c . For a con- 
tinuous (second-order) phase transition at r c , a quantum criti- 
cal point (QCP) emerges, which has wide reaching influence 
on the temperature dependence of thermodynamic and trans- 
port properties. In metallic systems near magnetic QCPs pro- 
nounced non-Fermi liquid (NFL) effects have been observed, 
which depend on the dimensionality (D), type of magnetic in- 
teraction, i.e. antiferromagnetic (AF) or ferromagnetic (FM), 
as well as type of quantum criticahty (itinerant, locally crit- 
ical, Kondo-breakdown, etc.) fl, S E B S i] • Such NFL 
effects are related to an anomalous enhancement of the en- 
tropy S(T,r) near the QCP, resulting from quantum critical 
fluctuations. For a pressure-tuned QCP (r oc p - p c , p c : crit- 
ical pressure) it has been pointed out, that the Griineisen ra- 
tio F oc a/C of thermal expansion a to specific heat C di- 
verges in the approach of the QCP II 1 OL 1 1 ill , as recently found 
for various systems Il2l Il3l Il4, llal . When the control pa- 
rameter to tune the system to the QCP is the magnetic field 
(r = H - H c , H c : critical magnetic field), a corresponding 
divergence has been predicted for the magnetic Griineisen pa- 
rameter T mag = -(dM/dT) /C 111 ill . This property equals the 
magnetocaloric effect, i.e is proportional to the slope of isen- 
tropes in the temperature vs. magnetic field phase diagram 
jl ill . Below, we demonstrate the use of magnetization for 
characterizing a NFL and report for the first time a divergence 
of Tmag in the approach of a (field-induced) QCP. 

We focus on tetragonal YbRli2Si2, which is a clean and 
stoichiometric heavy-fermion metal that displays a magnetic 
field-tuned QCP at H c = 0.06 T (0.66 T) for H _L c 
(H || c) lfl6ll . This QCP arises when very weak AF order- 
ing at Tn = 70 mK, with an ordered moment as small as 
2 x 10~ 3 jUft/Yb [17], is continuously suppressed by magnetic 



field. Various thermodynamic, magnetic and transport exper- 
iments on YbRli2Si2 as well as its slightly Ge-doped variant 
YbRh 2 (Sii-. v Ge A .)2 (T N = 20 mK, H c = 0.027 T, H ± c) have 
revealed evidence for quantum criticality, that is controlled 
by magnetic field JH, [H: At H - H c , C(T)/T exhibits 
a stronger than logarithmic divergence while p has T-linear 
temperature dependence. For H > H c , Fermi-liquid (FL) be- 
havior is induced, as evidenced from the electrical resistivity, 
described by p(T) - po + AT 2 , with the coefficient A(H) di- 
verging towards H c proving a field-induced QCP |16, 18]. The 
Sommerfeld coefficient y{H) in the FL regime (H > H c ) dis- 
plays a (H - H c y 1 ^ divergence, which is incompatible with 
the itinerant theory for a spin-density- wave QCP iflill . Bulk 
susceptibility [ 19, 20] and nuclear magnetic resonance exper- 
iments [21] indicate that in a wide regime of the T - H phase 
diagram the critical fluctuations have a FM character. Only 
for temperatures below 0.4 K and fields below 0.25 T, i.e. very 
close to the AF ordered phase, AF fluctuations dominate. For 
H > H c , a strongly enhanced Sommerfeld- Wilson ratio that 
exceeds a value of 30 upon approaching the critical field has 
been observed fl^l . Note, however, that the observed tem- 
perature dependence of the bulk susceptibility and spin-lattice 
relaxation rate in the quantum critical regime cannot be ex- 
plained within the itinerant theory for FM quantum critical 
fluctuations in either 2D or 3D [22]. 

The critical Griineisen ratio r„- = V mo i/Kr X a cr /C cr (V mo /: 
molar volume, kj\ isothermal compressibility), where a cr and 
C cr denote the volume thermal expansion and specific heat 
after subtraction of non-critical contributions BlOll has been 
studied for YbRli2(Sii_ A Ge A )2 at zero magnetic field and tem- 
peratures down to 80 mK [12]. Below 0.6 K, a T cr oc T' 01 
divergence has been obtained, which may indicate a local 
type of QCP IU2I1 . Indeed, detailed studies of the evolution 
of the Hall coefficient upon field-tuning the system suggest a 
drastic change of the Fermi volume due to a localization of 
the 4/-electrons at the QCP [23]. Furthermore, thermody- 
namic evidence for an additional energy scale T*(H), possi- 
bly related to the Kondo-breakdown was obtained from mag- 
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FIG. 1 : Magnetization divided by field M/H of YbRh 2 Si2 as a func- 
tion of temperature. Inset: Inverse of M/H vs temperature. Dashed 
line indicates linear fit for low-temperature region (T < 0.3 K). 



netostriction and magnetization experiments [22]. Below, we 
address the nature of quantum criticality and the evolution of 
the entropy across T*(H) by means of the magnetic Griineisen 
ratio, which is the appropriate thermodynamic property for a 
field-induced QCP. 

For our measurements, we used hi gh-q uality single crystals 
(po = 1 yuQcm) characterized before 11611 . The magnetization 
was measured utilizing a high-resolution capacitive Faraday 
magnetometer 12411 . Specific-heat measurements have been 
performed with the quasi-adiabatic heat-pulse technique. 

Figure 1 shows the magnetization divided by field M/H as 
a function of temperature. For fields above 0.1 T, M(T)/H 
tends to saturate, while it keeps increasing down to the lowest 
temperature below 0. 1 T. The strong increase of M(T) /H with 
decreasing temperature at the critical field // C =0.06T is re- 
flected in the negative curvature in the inverse of M/H (see the 
inset). The deviation from the Curie-Weiss law above 0.3 K 
may be due to FM correlations in this system. Linear fitting 
of the data in the low-temperature region (T < 0.3 K) yields a 
Weiss temperature 9 = -0.46 K and an effective magnetic mo- 
ment yU e ff=1.6jUB/Yt> 3+ ion, close to the reported values at zero 
field in a previous ac-susceptibility study itlol . At H > H c , 
M(T)/H saturates and does not reach zero. At a FM QCP, the 
magnetization is expected to diverge as T — >0 with 6-0. The 
finite negative 9 at H c indicates that AF fluctuations dominate 
close to the QCPjUl. 

The temperature derivative of the magnetization for vari- 
ous different magnetic fields at and above the critical field has 
been investigated between 70 mK and 3 K. In Figure 2, we 
compare the temperature dependence of -{dM/dT)/T, dis- 
played in part (a) with that of the respective specific heat co- 
efficient C(T)/T in part (b) and analyze the ratio between both 



quantities, T mag (Fig. 2c). At the critical field H c = 0.06 T, 
both quantities diverge upon cooling, but the divergence in 
-(dM/dT)/T is much stronger, resulting in a divergent mag- 
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FIG. 2: (Color online) (a) Temperature derivative of the magne- 
tization as -(dM/dT)/T (left axis), (b) electronic specific heat as 
C(T)/T and (c) magnetic Griineisen ratio r ma6 vs temperature (on a 
logarithmic scale) for YbRh 2 Si2 in magnetic fields applied perpen- 
dicular to c-axis. Volume thermal expansion of YbRh2(Sio.95Geo.o5)2 
as —fi/T at zero field is shown for comparison in (a) (right axis) 1 12]. 
Inset in (a): Expanded plot for -(dM/dT)/T vs T for #=0.3 T. The 
arrow indicates the position of the maximum. Inset in (b): r mag vs 
temperature for YbRh 2 Si2 on a log-log plot. The solid lines represent 
r ra a g (r) oc T e with 6=0.7 and 2.0 for low- and high-temperature re- 
gions, respectively. Inset in (c): Field dependence of saturated r raag , 
as T — > 0. The solid line indicates -G, {H - H{."'y l with G, = -0.3 
and H{" =0.065 T. 



netic Griineisen ratio. At H = 0.1 T, -(dM/dT)/T diverges 
even stronger than at 0.06 T in the investigated temperature 
regime. However, the analysis of the field-dependence of the 
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magnetic Gruneisen ratio, discussed below, is consistent with 
a QCP at 0.06 T and suggests a saturation of the 0.1 T data 
due to a crossover from NFL to FL behavior at temperatures 
below 70 mK. Such crossover scales are known to be different 
for different physical quantities [5]. Indeed, the specific heat 
coefficient C(T)/T displays such a NFL to FL crossover for 
H = 0.1 T already around 100 mK, whereas the respective 
crossover in electrical resistivity at the same field is located at 

80 mK Cl- 
in Fig. 2(a) (right axis), we also show the volume ther- 
mal expansion as /3/T for YbRli2(Sio.95Geo.o5)2 01211 for com- 
parison. The divergence of -(dM/dT)/T is stronger at 
high temperatures and becomes comparable to /3/T at low 
temperatures. The ratio between f3 and dM/dT gives the 
pressure dependence of the critical field as expressed by 
dH c /dp=V„,/3/(dM/dT), where V m is the molar volume ill ill . 
This ratio is temperature dependent and decreases with de- 
creasing temperature (not shown). At low temperatures it 
tends to saturate and the extrapolated value of V m /3/ (dM/dT) 
for T -» is (0.19 + 0.05)T/GPa which is consistent with 
previous hydrostatic pressure experiments lf20ll . 

Interestingly, -(dM/dT) /T for H > 0.2 T shows a max- 
imum (cf. inset of Fig. 2(a) for data at 0.3 T) whose posi- 
tion coincides well with the temperature where the specific 
heat coefficient passes a maximum, T max . Similar behavior in 
C(T)/T has also been observed in other heavy-fermion com- 
pounds such as CeCu5 9Auo i and Ceo.sYo.2Cu2Si2 1I25I l26ll . 
The origin of such maxima may be understood qualitatively 
by the Zeeman splitting of the Kondo resonance 1 27] which 
results in a more polarized state. Our results as well as those 
for CeCu 5 . 9 Auo.i and Ceo.8Yo.2Cu 2 Si2 10 IH show that the 
maximum temperature in C(T)/T increases linearly with in- 
creasing field. For the magnetic Gruneisen ratio F mag (r), the 
maxima in C/T and -(dM/dT)/T cancel each other, leading 
to a monotonous crossover from NFL to FL behavior at fields 
larger than 0.2 T (Fig. 2(c)). 

Before discussing details of F mag (r) data, we first summa- 
rize the important theoretical conclusions most related to the 
present study with respect to the degree of assumptions. 

(1) r mag diverges in the approach of any field-induced QCP 
whenever the characteristic energy scale of a system is con- 
tinuously suppressed to absolute zero fioll . 

(2) Assuming scaling, the critical behavior is governed by 
the correlation length £ and the correlation time £ T such that 
the temperature dependence in the quantum critical regime 
T m!lg (T,H = H c ) oc T- 1/vz (v and z are the correlation-length 
exponent, and dynamical critical exponent, respectively). For 
the field dependence in the FL regime, scaling analysis re- 
markably predicts not only a universal functional dependence 
but even its prefactor without any adjustable parameter, i.e. 
T lmg (T = 0,H)=-G r (H - H c )-\ with G r = v(d- z) (d: di- 
mensionality of the critical fluctuations) HTot [Till . Further- 
more, this prefactor equals the exponent in the divergence of 
the Sommerfeld coefficient 7 oc (H-H C ) G ' JH[iT|]. Thus, the 
character of the QCP is completely determined by the values 
of v, z and d. 
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FIG. 3: (Color online) Magnetization difference divided by tempera- 
ture increment -AM/ AT = -{M(T + AT,H) - M(T - AT,H)\/2AT 
vs magnetic field for YbRh2Si2, obtained using isothermal magneti- 
zation measurements at T = 0.08 K (black circles), 0.33 K (red open 
squares), 0.75 K (blue open circles) and 1.5 K (green squares). Inset: 
H-T phase diagram of YbRh 2 S12 with the peak positions in - AM/ AT 
(black solid circles) and T*(H) (grey solid line) as determined from 
transport and thermodynamic properties (red symbols) 1 22, 23|]. The 
black and blue lines represent the phase boundary of the antiferro- 
magnetic (AFM) ground state and crossover to the Fermi liquid (FL) 
regime, respectively 12211 . 



(3) Within the itinerant theory, the correlation length exponent 
v - 1/2 and the dynamical critical exponent z equals 2 and 3 
for AF and FM case, respectively Jj]. 

We now check the consistency of our experimental results 
with these theoretical predictions. Obviously F maa diverges 
as a function of temperature at H-H C -Q.Q6T (Fig. 2(c)) 
proving the existence of a field-induced QCP (assumption 
1). The saturation values of r mag for T — > as a func- 
tion of field, plotted in the inset of Fig. 2(c), diverge like 
(H - H c y l as expected from scaling (assumption 2). The 
fit reveals -G r (H - fff') -1 with fff'=(0.065+0.004)T and 
G r = -0.30 + 0.01. For YbRh 2 (Sio.95Ge o .05)2, the field de- 
pendence of y within the FL regime has been studied in detail 
revealing y oc (H - H c )~ 033 H 1 8fl which has an exponent very 
close to our value of G r . This underlines the thermodynamic 
consistency of the data. We note also that the obtained value 
of G r sets strong constraints for the scaling parameters v, d 
and z. Within the itinerant theory (assumption 3), in which 
v - 1/2, the observed G r would require d — z = 2/3 which is 
impossible. By contrast, a critical Fermi surface model which 
may be relevant for a Kondo-breakdown QCP has been pro- 
posed |28J], in which the electronic criticality is described by 
v = 2/3, z = 3/2, and d = 1, yielding G, = -1/3 similar is in 
our experiments. 

Next, we focus on the temperature dependence of the mag- 
netic Gruneisen ratio within the quantum critical regime. In- 
terestingly, divergent behavior is found not only at H c but 
also at 0.1 and 0.15 T (see Fig. 2c). However, the obtained 
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flf '=0.065 T is very close to 0.06 T, supporting the QCP at 
H C -Q.Q6 T. Furthermore, the extrapolated values of the sat- 
urated r mag at ff=0.1T and 0.15T from -G r (H - fff')" 1 
are 11.8 and 26.6T -1 , respectively, strongly exceeding the 
largest measured values shown in Fig. 2(c). Thus, the low- 
est temperature in our study (~70mK) is not low enough to 
observe saturation of r ma£ , at fields very close to the QCP and 
thereby explains the absence of saturation in -(dM/dT)/T at 
0.1 and 0.15 T. As shown in the double-log plot in the inset 
of Fig. 2(b), r mag oc T e with e « 2.0 from 3 K down to 
0.3 K. This temperature dependence is much stronger than 
any expectation within the itinerant theory. At lower tem- 
peratures, a crossover to a weaker divergence is found with 
e « 0.7 below 0.25 K, which interestingly is similar as the 
exponent obtained by the study of the thermal Griineisen ra- 
tio for YbRli2(Sio.95Geo.o5)2 at zero field [12]. The crossover 
around 0.3 K may be related to the interplay of AF and FM 
fluctuations, as it is of the same order of magnitude as the ab- 
solute value of the Weiss temperature lfl9ll . The fact that no 
single power-law behavior is found over a larger temperature 
regime makes the interpretation of the Griineisen ratio within 
the quantum critical regime difficult. 

At last, we study the evolution of the magnetic entropy 
upon crossing T*(H). For this purpose, it is important to 
follow the field dependence of dM/dT, which according to 
the Maxwell relation equals the field dependence of dS / dH. 
Using isothermal magnetization scans at different tempera- 
tures, we calculate the difference in magnetization divided by 
the temperature increment, -AM /AT = -{M(T + AT, H) - 
M(T - AT, H)}/2AT vs magnetic field, plotted in Fig. 3. For 
each curve, we have used two sets of isothermal magnetiza- 
tion M(H) data at temperatures T - AT and T + AT. The 
-AM / AT(H) traces show a peak whose position shifts to 
larger field with increasing temperature as displayed in the 
temperature versus field phase diagram in the inset of Fig. 3. 
The peak positions agree satisfactorily with the energy scale 
T*{H) found previously in the Hall coefficient fl23ll . magne- 
tostriction and susceptibility 02211 (cf. red symbols in the inset 
of Fig. 3). Through AM/ AT * dM/dT = dS /dH, the field de- 
pendence of the entropy S (H) could be obtained by integrat- 
ing AM I AT over the field. Since dM/dT < in the entire field 
range (at all accessible temperatures), the magnetic entropy 
decreases with increasing field. The positions of the maxima 
in -AM/ AT correspond to inflection points in S(H), i.e. they 
mark characteristic magnetic fields, at which entropy S (H) at 
constant temperature is reduced most strongly with increasing 
field. Remarkably, these inflection points of the entropy agree 
very well with T*{H) determined from Hall effect 12311 . mag- 
netostriction and magnetization 12211 and thus provide further 
thermodynamic confirmation of this additional energy scale in 
YbRh 2 Si 2 . 

In conclusion, our study of the low-temperature magne- 
tization and specific heat of YbRli2Si2 reveals for the first 
time in any system a divergence of the magnetic Griineisen 
ratio F mag = -(dM/dT)/C at a magnetic field-tuned QCP. 



This property provides information on the scaling of the mag- 
netic entropy due to quantum critical fluctuations. The field- 
dependence of the magnetic Griineisen ratio within the FL 
regime follows F mag (r -» 0) = -G r (H - H c y l [0 with a 
prefactor G r = -0.3 + 0.01 that is consistent with a recent 



critical Fermi surface model 12811 . Within the quantum critical 
regime, complicated behavior with a crossover scale around 
0.3 K is found which may be related to the interplay of AF 
and FM fluctuations in this system. We also observe an in- 
flection point in the entropy at the scale T*(H) previously ob- 
served in Hall- and thermodynamic measurements 11221 12311 . A 
reduction of spin entropy when entering the entangled heavy- 
electron fluid upon increasing the magnetic field is consis- 
tent with the derealization of /-electrons upon field tuning 
through the QCP. 
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